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Abstract. Clifton's exact solution of f{R) = iJ^^* gravity describing 
a dynamical spherical metric which is asymptotically Friedmann-Lemaitre- 
Robertson- Walker is studied. It is shown that it harbours a strong spacetime 
singularity at a finite radius and that this singularity is naked. 
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1. Introduction 

Type la supernovae have provided us with the knowledge that the universe is currently 
in a phase of accelerated expansion yy . This acceleration has been modelled in various 
ways; the most common models are probably dark energy ones based on General 
Relativity (hereafter GR, see f5] for a list of references). However, the exotic and ad 
hoc dark energy leaves many cosmologists dissatisfied and attempts have been made 
to model the cosmic acceleration without dark energy. f{R) theories of gravity akin to 
the quadratic theories required by the renormalization of GR have been introduced in 
the metric 3 , Palatini and metric-affine ^ formulations and have received much 
attention in recent years (see [B] for a review and [7] for short introductions). 

Along with cosmological and other considerations {e.g., stability, weak-field limit, 
ghost content), it is important to understand spherically symmetric solutions in these 
theories, a task which has proved to be non-trivial (see [8] and references therein). For 
definiteness, we consider here metric f{R) gravity described by the action 

S=^ [ d'^xy/^fiR) + S""""'^'^ , (1) 
2k J 

where f{R) is a non-linear function of its argument and 5'™"**'^'' is the matter part of 
the action. R denotes the Ricci scalar of the metric gab with determinant g, k — SttG 
where G is Newton's constant, and we adopt the notations of Ref. ^. 

It is well known that the Jebsen-Birkhoff theorem does not hold in these theories, 
which adds to the richness and variety of spherically symmetric solutions. Of particular 
interest are black holes in generalized gravity, which have been studied especially in 
relation to their thermodynamics Q {e.g., [ll')- Since f{R) theories are designed 

I The thermodynamics of local Rindler horizons in f{R) gravity, which is used to derive the classical 
field equations as an equation of state [TO) is modeled after the thermodynamics of dynamical f{R) 
black holes. 
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to produce a time-varying effective cosmological constant, the spherically symmetric 
and black hole solutions of interest likely represent central objects embedded in 
cosmological backgrounds. Not much is known about this kind of objects even in the 
context of GR, although a few examples are available P ^[rai[TIi ni l [ro i [T7 1 [T5 1. [ra] 
in Einstein's theory. Even less is known about f{R) black holes and spherically 
symmetric solutions, which deserve to be understood better. Here we consider 
a specific solution proposed in /(i?) = _R^+* gravity in [20]. The observational 
constraints set the limits 6 = (— 1.1± 1.2) • 10^^ on the parameter 5 [20l[2T], while 
local stability requires f"{R) > [551 [23j, i-c, S > hence we will use positive values 
of this parameter. 

The solution proposed in [20] is dynamical and presumably represents some kind 
of dynamical central object embedded in a spatially flat FLRW background in vacuum 
f{R) = R^~^^ gravity. This solution is made possible by the fact that the fourth order 
field equations of vacuum metric f{R) gravity 

f'{R)Rab -^9ab^ V,V,/'(i?) - gatafiR) (2) 

can be rewritten in the form of effective Einstein equations with geometric terms acting 
as a form of effective matter as 



1 1 

Rab — 77 9abR 



^a^bf - gabOf + gab ^ 



(3) 

In this picture the effective matter spoils the Jebsen-Birkhoff theorem and fuels the 
cosmic acceleration. Alternatively, an equivalent representation of f{R) gravity as a 
Brans-Dicke theory with a scalar field potential exhibits a massive spin zero degree 
of freedom that causes these effects ^6j. Since exact spherically symmetric dynamical 
solutions of f{R) gravity in asymptotically FLRW backgrounds are harder to find than 
in GR (where few are known anyway) and are therefore valuable, we study Clifton's 
solution in the following. 

2. Clifton's spherically symmetric dynamical solution 

Clifton's spherically symmetric dynamical solution in vacuum fiR) = R^+^ gravity 
[20] is given by 

ds^ = -A2ir)dt^ + a^{t)B2{r) {dr^ + r^dn^) , (4) 
where d^^ = d6^ + sin^ 6 dip"^ is the line element on the unit 2-sphere, 

B2{r) = + A2{r)'^+^'-^ , (6) 



r 

i5(l + 2i5) 

a{t) = t^^ , (7) 

q'^ = 1 - 2(5 + 4(5^ , (8) 

in isotropic coordinates and using the notation of [20 for the metric components. The 
line element (|4]) reduces to the Friedmann-Lemaitre- Robertson- Walker one in the limit 
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C2 ^ 0. In the limit (5 — > in which the theory reduces to GR, the metric (H]) reduces 
to the Schwarzschild solution in isotropic coordinates. This suggests that the positive 
root be taken in the expression q = ±Vl — 2S + 4(5^ deriving from eq. ([8]), and that 
g ~ 1 — (5 in the limit \6\ << 1. Moreover, only positive values of the constant C2 will 
be considered since the latter reduces to the Schwarzschild mass in the limit to GR. 
The solution ([l])-!!]) is conformal to the Fonarev solution [24] which is conformally 
static |25j . and therefore is also conformally static. This is a property shared with the 
Sultana-Dyer solution [12| and with some representatives of the class of generalized 
McVittie solutions ^ . 

We now want to write the metric (|4]) in the Nolan gauge, in which it is 
straightforward to identify the apparent horizons (if they exist). To this end, we 
make use of the Schwarzschild-like radial coordinate 

1 + — I , (9) 



C2 
r 



in terms of which dr = ( 1 — 



df and we eventually transform to the areal radius 



R: 



a{t)^B2ir)r 



a{t)fA2{rY 



The line element (jlj then becomes 

ds^ = -A2dt^ + a^Af-^df^ + R^dVt^ 
Using the fact that 



df 



dR-A^ 



afdt 



dR-A^ " afdt 



\^ 2 _|_ 2(g+2<5-l) Ci /i^ 2 



-C[r) 



(10) 



(11) 



(12) 



where an overdot denotes differentiation with respect to t and 



c(0 = i+^^^±^^A-^i- 

q f 

the metric assumes the Painleve-GuUstrand-like form 

q + 2(25-l) 

A 



2(g+2(5- 1) C^a = 



A. 



1 



A, 



AiC^ 



■ 2 ~2 

■ a r 



dt^ - 



2A 



AiC^ 



-j^A ' (13) 



'•dtdR 



dR^ 



R^dn'' 



In order to eliminate the cross-term in dtdR we introduce the new time coordinate t 
defined by 

1 



dt = 



[dt + Pit, R)dR] , 



F{t,R) 

where F(t, R) is an integrating factor that satisfies the equation 
dR\F dt\F, 



(15) 



(16) 



Clifton's spherical solution in f{R) vacuo harbours a naked singularity 



to ensure that dt is an exact differential. The line element then becomes 

q + 2(2S-l) 



ds^ = - An 



1 - 



AiC^ 



■ 1 -2 

■ a r 



F^dP 



+ 2F I A2(3 



q+2(25-l) 



A, 



AiC^ 



■ 2 ~2 

■ a r 



A. 



A'iC'^ 



df > dtdR 



+ 



A, 



AiC^ 



■ 2 -2 

■ a r 



^2 ^ 2^2 



AiC^ 



■ dfP ■ 



AiC^ 



By setting 



q+2S-l 



ar 



AlC 



<;+2(25-l) 
\2 



■ cfir^ 



(17) 



(18) 



the dtdR cross-term disappears and we arc left with the Nolan gauge metric 

2q(\-q)J^2S(2-q)-2 

1 



ds2 ^ _ A2DF'^df 
+ R^dfl'^ , 



AlC^ 



1 + 



dR'' 



where H = d/a and 



9+2(2^-1) 

A " 
D=l ^2 



g(2-q) + 2a(2-<,)-2 



AlC^ 



■ 2-2 1 

a r = 1 — 



AlC^ 



H^R^ . 



(19) 



(20) 



The apparent horizons, if they exist, are located at = 0. This equation is satisfied 
if A2 = 0, or £) = 0, or C = 0, corresponding to 



r = C2 



q(2-q)+2S{2-q)-2 

AlC^ = A2 " H^R'' , 



(21) 
(22) 



2{q + 2S-l) C2A 2^ 







(23) 



respectively. The locus r = C2 for which A2 — (which describes the Schwarzschild 
horizon in the limit 5 ^ in which the theory reduces to GR) corresponds to a 
spacetime singularity. In fact, the Ricci scalar is 



DO 

a 



A2{r) 



(24) 
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and diverges as r ^ C2 (it reduces to the familiar value 6 + 2H j in the 

0-2^0 limit). Furthermore, this singularity is a strong one in the sense of Tipler's 
classification 26 : the metric determinant is 



g = -a''(tyi^l + ^j A2{ry^+'''-' (25) 

and vanishes as r — > C2. The volume of a body is shrunk to zero as it approaches 
this singularity and the energy density of a (real or effective) fluid diverges there. No 
object can cross the locus r = C2. The regions < r < C2 and r > C2 describe two 
disconnected spacetimes. It seems that the pull of the effective matter in the universe 
has stretched the r = singularity of the Schwarzschild black hole into a sphere. [§| 

Eq. (pS]) . corresponding to C = 0, has no solutions for 5 > 0, C2 > 0, and R> 0. 
In fact, for < (5 << 1, this reduces to 1 + 2(5 '^g. ^ = 0, which cannot be satisfied 

r{l--^) 

if C2rS > 0. Let us focus on eq. (|22|) corresponding to £> = 0. This yields 



HR = ± 



2(<z + 2<5- 1) C2a 



q R 



^2 



A, ' . (26) 



In the limit of small S this equation reduces to HR — ± 



1 I 2SC2a /i"(l"¥) 



and, in an expanding universe in which HR > 0, we discard the negative sign in 
eq. (pS)) . The apparent horizons, if they exist, are located at the roots of the equation 

2g^ -2q-45 + 2(5q + 2 Ofrr-UOX 1\ 2g^ - 2q -45 + 2(5g + 2 , 25-1-g 

HR^ - A, ~- R - ^^'^^^^ C2aA, ~^ 

q 

= . (27) 

Although eq. (j27p is written in the form of a quadratic algebraic equation, it is really an 
implicit equation for the i?-coordinate of the apparent horizons because the coefficients 
are functions of r{R). In spite of this fact, it is still useful to regard eq. (|27|) as a formal 
algebraic equation. 

To gain some insight, consider the following two limits. In the limit C2 — > in 
which the central object disappears and the solution is a FLRW space, r — r becomes 
a comoving radius and R becomes a proper radius, while eq. (|27p reduces to 

RiHR-l) = (28) 

which yields as a solution Rc = 1/H, the radius of the cosmological horizon. 
In the limit (5 — s- in which the theory reduces to GR, the exponent 

2q^ -2q-4:5 + 2Sq + 2 



2q 



1 - (5 ^ 1 , (29) 



while -^(g+^'^-i) K -26 ^ 0, H ^ and eq. ^ simply yields A2R = 0, or 

afA2 ^ =0 and A2 = 0, hence r — C2- This is the usual Schwarzschild 
horizon expressed using the isotropic radius (the corresponding Schwarzschild radius 
is f = 2C2). 



§ This feature is consistent with the known phenomenology of the Sultana-Dyer solution, of 
generalized McVittie solutions |27| , and of higher-dimensional Gauss-Bonnet black holes |17| . 
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Returning to the general case, we see that eq. (^7)) has the formal solutions 

2,2 



^2 

Ri.2 — — 



-2g-4.S + 25g + 2 



2H 



2q'^ -2q~45 + 25q + 2 

2H 



.(30) 



The root Ri corresponding to the positive sign in eq. ([50)1 yields a cosmological horizon. 
For the physical range of parameters < 5 << 1 and C2 > 0, in an expanding 
universe (a > 0) the argument of the square root is larger than unity and the root R2 
corresponding to the negative sign in eq. (j30[) is negative and unphysical. Therefore, 
we conclude that there is no black hole apparent horizon and the singularity at r = C2 
(or i? = 0) is naked. 

For completeness, we can consider also a contracting universe with a < 0, which is 
obtained for (5 < 0, although this situation is clearly not interesting for f{R) theories 
aiming at explaining the current acceleration of the universe and <5 > is required to 
stabilize the theory against explosive local instabilities [23]. In this case the negative 
sign has to be chosen in eq. (P^ . leading to 



HR^ 



A, 



-2q-45 + 2(5q + 2 

R 



2{q + 25-l) 



Ay 



= . 

The formal solutions are 



R: 



3,4 



2q2 -2q-AS + 2Sq + 2 
A2 

2\H\ 



1± Wl 



(7 + 2(5-1 



C2 |a| a: 



2i5-l-g 



(31) 



(32) 



and are both non-negative. When the argument of the square root is positive the 
upper sign yields again a cosmological horizon while the lower sign yields a black hole 
apparent horizon. In the solution (HI)-® the scale factor a{t) ~ l/i'"'' has a pole-like 
singularity at i = and a is always negative, ensuring the existence of the black 
hole apparent horizon at all times. However, as already remarked, this situation is 
completely unphysical. 

3. Quasi-local mass 



The mass of the naked singularity is also of some interest. Due to the fact that 
the solution (H])-® is not asymptotically flat, the ADM mass is not defined and one 
needs to resort to the concept of quasi-local energy on a 2-surface surrounding the 
singularity. Thanks to the spherical symmetry of this solution it is straightforward to 
compute the Hawking-Hayward [28l [29] and the Misner-Sharp [30] energies. 
By introducing the affine parameters ^ and 77 according to 



yJMDFdt- 



^/A^Fdt 



1 



Ate 



1 



A^H^R^ 



1 + 



C^D 



dR 



(33) 



(34) 
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^^Ml-,) + 2Si2-,)-2^ (35) 



with 



the hne element ()19p is rewritten in the standard form 

ds'^ = -2d^dr] + R^dn^ (36) 

and the Hawking-Hayward quasi-local energy Mhh on a 2-surface of constant radius 
R is given by [29| 

M„„ . R iR,R, + i) ^ f (l - c^ff^H^H^ ^ m 

In the presence of spherical symmetry the Misncr-Sharp quasi-local energy Mms 
on a 2-sphere R =constant is defined by [30] 

l-^^=-V=i?V,i?, (38) 
it 



which yields 

2 V C^D + APH^R^^ 

In the case of a contracting universe, the two mass notions coincide on the black hole 
apparent horizon Rah given by £> = 0, i.e., Mhh = Mms = Rah/2. 



4. Conclusions 



In view of the fact that cosmology may be showing us the first-ever detected 
deviations from Einstein's gravity and of the attention given to f{R) gravity theories as 
possible models of the cosmic acceleration, it is of great interest to understand black 
holes and other spherically symmetric solutions of f{R) gravity. Since the Jebsen- 
Birkhoff theorem does not hold in these theories, spherically symmetric solutions 
do not have to be static. These theories are designed to produce an eflFective 
dynamical cosmological constant to reproduce the current acceleration of the universe 
and, therefore, dynamical exact solutions with spherical symmetry describing a 
central object embedded in a cosmological background are particularly valuable. 
Unfortunately, this kind of solutions is poorly understood already in the context 
of GR and deserves more attention in the future. Few examples are available 
[Tl[T31[Tl[Tni[Tni[Il[ISl[IS] and further work is in progress on generalized McVittie 
solutions [271. In particular, it seems difficult to find generic black holes solutions 
embedded in cosmological backgrounds. The goal of finding interior solutions for 
spherically symmetric f{R) gravity (mainly with numerical methods) seems also a 
worthy one [8]. All these issues will be addressed in future publications. 
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